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ABSTRACT 

We report the discovery of a pair of infrared, axisymmetric rings in the planetary 
nebula NGC 1514 during the course of the WISE all-sky mid-infrared survey. Simi- 
lar structures are seen at visible wavelengths in objects such as the "Engraved Hour- 
glass Nebula" (MyCn 18) and the "Southern Crab Nebula" (Hen 2-104). However, in 
NGC 1514 we see only a single pair of rings and they are easily observed only in the 
mid-infrared. These rings are roughly 0.2 pc in diameter, are separated by 0.05 pc, 
and are dominated by dust emission with a characteristic temperature of 160 K. We 
compare the morphology and color of the rings to the other nebular structures seen at 
visible, far-infrared, and radio wavelengths, and close with a discussion of a physical 
model and formation scenario for NGC 1514. 

Subject headings: infrared:stars, planetary nebulae: individual (NGC 1514), surveys:WISE 
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1. Introduction 

In 1790, William Herschel found "a most sin- 
gular phaenomenon! A star of about the 8th 
magnitude, with a faint luminous atmosphere, 
of a circular form, and of about 3' in diameter." 
This object, now known to us as NGC 1514, 
with its central star inseparable from a "shining 
fluid" , convinced him that not all the nebulae he 
studied could be r esolved into clu sters of stars 
as he had thought ( Herschellll79ll ). Thus began 
over 200 years of study of NGC 1514, though it 
is still classified as the round or elliptical plan- 
etary nebula that Herschel saw. 

In modern terms, NGC 1514 is a moder- 
ately high excitation planetary nebula (PN) 
in Taurus (a = 04 h 09 m 16*990, 5 = +30° 
46' 33"44, J2000). The optically visible cen- 
tral source (CSPN) of NGC 1514 is unusual in 
that its spectral type is too late (AO) t o pro - 



duc e the observed nebu la. Kohoutek ( 19671 ) 



and iKohoutek & Hekelal (jl967n ~ developed the 



notion of a binary central source, and concluded 
that an My = —1.4 AOIII giant and a sub- 
luminous, dwarf O star would repro duce the ob 



served UV spectrum. By contrast, iGreenstein 



concluded through the use of higher 
spectral resolution spectrophotometry that the 
visible star is a horizontal branch A star. This 
star is slightly cooler and has a much lower 
absolute magnitude than an AOIII giant. He 
finds that the CSPN is better represented by 
an My = +0.8 A star of < 10000 K, and an 
sdO star with a tempe rature of 100,000 K and 



My = +2.8. Later, iFeibelmanl (119971 ) built 



on earlier IUE observations by Seatonl (jl980l ) 
to show that the CSPN could be represented 
by an early A star of ~ 9000 K along with a 
> 60,000 K subdwarf; he also found that the 
CSPN's ultraviolet flux varied by more than a 
factor of two on timescales of a year to a decade. 

NGC 1514 itself is typically catalogued as 
a round or slightly elliptical nebula , with 



an amorphous appearance (e.g. iBalickl 119871 ). 
Kohoutek! (jl968l ) describes the nebula as having 
an inner shell (or main body) with a diameter of 
~ 136", and an outer, homogeneous, spherical 



layer of ~ 206". He interprets condensations 
within the inner part of the nebula toroid 
with t he axis of symme try at a position angle 
of 35°. IChu et al.1 (119871 ) classify it as a double 
shell PN, 132" and 183" in diameter. Figure Q] 
illustrates the optical appearance and uses their 
nomenclature (which we adopt) to highlight the 
various features. 
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Fig. 1. — Labeled, logarithmically-scaled, Ha 
image of NGC 1514. (FITS data kindly supplied 
by R. Corradi.) 



Hajian et al 



(119971 ) describe NGC 1514 as 
a "lumpy nebula composed of numerous small 
bubbles" based on their deep [O III] image. 
They interpret the bubbles at the edge of the in- 
ner shell to be sweeping up outer shell ("halo") 
material which is providing pressure confine- 
ment. They confirm a somewhat filamentary 
structure in the outer shell t hat h ad been re- 
ported earlier by IChu et al.l (|1987l ). Unfortu- 
nately, NGC 1514 does not have an observable 
halo of the type needed for the timescale corre- 
lation method for determining the distance that 
is the focus of their paper. 

Per haps the spatiokinematic study by Muthu &; Anandaraol 
(j2003l ) is the most in-depth look at the forma- 
tion and structure of NGC 1514. They used an 
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imaging Fabry-Perot spectrometer to produce 
velocity-resolved maps of the 5007 A [O III] 
line. From modeling the double- and triple- 
component line profiles seen in these maps, 
they deduced that NGC 1514 is a generally 
ellipsoidal shell nebula, but that it also has two 
"blobs" or bubbles of material emanating from 
the center. The southeast bubble is slightly 
blue-shifted, and the northwest blob is red- 
shifted, thus defining a "polar" axis. The bub- 
bles themselves are uncollimated in that there 
are large velocity dispersions across each of the 
bubbles; this is in contrast to other highly col- 
limated, bipolar flows seen in many other PNe. 
Finally, they model the overall velocity struc- 
ture of the nebula and derive peak expansion 
velocity of 23 km s _1 near the center, falling off 
with radial distance as expected from projection 
effects along the line of sight. 

2. Observations 

The Wide-field Infrared Survey Explorer 
(WISE, IWright et al.1 l2O10h is a NASA Mid- 
class Explorer mission that is surveying the 
entire sky at wavelengths of 3.4, 4.6, 12, and 
22 /im (Wl through W4, respectively). The 
point source sensitivity (defined ciS ct 5<7 detec- 
tion with eight repeat exposures per source) 
is currently estimated to be 0.08, 0.11, 1, and 
6 mJy in these passbands, respectively. WISE 
surveys the sky in 47' wide strips with a signif- 
icant overlap between each strip that yields the 
eight or more coverages (depending on ecliptic 
latitude). NGC 1514 was thus observed twelve 
independent times on 2010 Feb 18-19 UT for 
a total of 106 seconds of integration time in 
each of the four bands; NGC 1514 was easily 
detected in all four bands. 

The data presented here were processed with 
the "first-pass" version of the WISE pipeline 
that utilized preliminary reduction algorithms 
and calibrations derived from early survey data. 
The twelve individual fields have been corrected 
for detector artifacts, optical distortions, and 
latent images. They were then combined into 
"ops coadds" with pixel outlier rejection to sup- 



press cosmic rays and transient noise events, 
and have been resampled from the native plate 
scales of 2.75 (Wl, W2, and W3) and 5.5 (W4) 
arcsec pixel" 1 to the WISE image atlas scale 
of 1.375 arcsec pixel -1 . The photometry has 
been calibrated using WISE measurements of 
blue (stellar) calibrators; calibration using red 
(ULIRG) calibrators currently differs by several 
percent, so that the absol ute fluxes are curren tly 
uncertain by ±10%. (See lWright et al.ll2010l for 
further discussion.) No color corrections have 
been applied to the photometry reported here. 
For our purposes, the correction terms for stel- 
lar spectra and blackbody continua and for all 
sources in the W4 band are essentially zero. In 
W3, the corrections do not exceed 15% until 
objects are cooler than ~ 160 K. For W2, the 
equivalent temperature is ~ 250 K, and for Wl, 
~ 400 K. 

The W4 data were also reprocessed using the 
deconvolution capabilities of AWAIC, a coaddi- 
tion and deconvolution package developed for 
WISE, but which can be used with generi c 
astronomical data sets ( Masci Sz Fowler 2009). 
The deconvolution algorithm "HiRes" is based 
on the Maximum Correlation Method, which 
was developed to enhance the scientific return 
from the IRAS survey. W4 data covering a 
24 arcminute field was passed through just three 
iterations of HiRes, enough so that the model 
image has an apparent spatial resolution equiv- 
alent to the native W3 resolution. The result- 
ing image, denoted as W4d eC on> is used only 
in morphological discussions; all flux measure- 
ments are based on the W4 ops coadd image. 

3. Results 

Inspection of even the raw, unprocessed 
WISE images immediately revealed an unex- 
pected structure in NGC 1514. A pair of bright 
axisymmetric rings that surround the visible 
nebula are plainly apparent (coadded data are 
shown in Figure [2^,); in fact, they dominate 
the emission at W4, giving the red appearance 
to the rings in the figure. Such a structure is 
not suggested in any of the visible wavelength 
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images available to us; e.g. the DPOSS data 
presented in Figure [21 the Ha image i n Fig - 
ure [p the p u blishe d figu res of iBalickl (11987). 
Haiian et all (Il997h . and ICorradi et al.1 (|2003l . 
and private communication), or the deep, very 
narrow-band image by D. GoldmarQ. 

Remarkably, the inner shell has almost ex- 
actly the same appearance in the mid-infrared 
as in the visible when the differing spatial reso- 
lutions are taken into account. Panels c and d 
of Figure [2] show that the inner shell is "gray" 
across all wavelengths (with the exception of the 
rings). The spatial extent, the positioning of 
the bubbles, the relative brightnesses, etc. are 
all identical at the ~ 10-20% level over many 
different wavelengths. This is in marked con- 
trast to other PNe such as those illustrated in 



Balick fc Frank! (|2002T ) and iHora et al.1 (12004). 
though the Group 1 PNe of IChu et al.1 (|2009l ) 
have 24 /im morphologies that compare favor- 
ably to the optical red or Ha images. Images of 
these objects at other wavelengths are needed 
to determine if they are similar to NGC 1514 in 
this regard. 

3.1. Photometry 

The spatial resolution of WISE does not al- 
low much new insight directly into the nature 
of the central source of NGC 1514; the 6" res- 
olution has no hope of resolving the central bi- 
nary. If the sdO+AOIII (or horizontal branch 
AO) binary model is correct, then at best we 
would expect the infrared fluxes to approxi- 
mately follow the Rayleigh-Jeans law for a ~ 
10,000 K blackbody. To verify this, we present 
infrared fluxes from the 2MASS and WISE sur- 
veys in Table [TJ These values are obtained 
by profile fit photometry as returned by the 
2MASS and WISE data pipe l ines. Encour- 
agingly, iTaranova &: Shenavrinl (|2007l ) got very 
similar magnitudes in the near-infrared using an 
aperture photometer: J = 8.238, H = 8.108, K 
= 8.036, and L = 7.927. 

The visual extinction to NGC 1514 has been 
estimated by various authors to be 1.5 to 2.5 



magnitud es, based primari ly on the E(B — V) 
color (e.g. Greensteinl 1972) or from the H a/H/3 
or radio/H/3 ratios (e.g. Cahn et all 19921 ). We 
use the new WISE fluxes to provide a much 
longer baseline in an effort to better estimate 
the correct visual magnitude of NGC 1514 pho- 
tometrically. We begin by constructing a model 
CSPN spectrum using visible-wavelength spec- 
tral templates for an Q5V star and an AOII I 
star (from STELIB, iLe Borgne et al.l 120031 ). 
and a 9800 K blackbody to represent the in- 
frared. The precise spectral types are relatively 
unimportant; we are primarily interested in 
correctly estimating the overall shape of the 
SED, not in the particulars of any spectral fea- 
tu res. Using the inters tellar extinction curves 
of iDraine &: Led (|1984l ). we redden the model 
until an acceptable match is achieved. (Extinc- 
tion curves of other authors were checked, but 
the differences are not large enou gh to affect 
our conclusions, and lDraine Led are the only 
set comprehensive enough to cover ultraviolet 
through infrared wavelengths.) 

The best model fit is shown in Figure (3j an 
extinction value of Ay = 1.6 is found. The fit is 
quite good except for the J band near-infrared 
point and the W3 mid-infrared point. Given 
the good agreement between the 2MASS pro- 
file fit value and the apertur e-derived value of 
Taranova &z Shenavrin ( 20071 ) . the 



' excess at 

J is real, at least in the sense that it is not a 
measurement error and it cannot be fit with 



http : //dg- imaging. astrodon. com/gallery/display . cf m? 



our si mplistic photospheric model. IWhitelock 
(|l985l ) forecast the great strength of the He I 
triplet at 1.083 /xm in m any PNe, while P a/? can 
also be very strong (e.g. Hora et al. 19991 ). so it 
is possible that some nebular emission is being 
picked up where it overlaps the CSPN. (Though 
to be fair, the K band will pick up Br7 and may 
also be affected.) The small excess at 12 fim is 
also real, and is definitely due to the nebular 
line emission (see Section [3.4p . 

Using this value of Ay , we deredden the mag- 
nitudes for the CSPN, and these are also re- 
ported in Table [TJ If we assume that the near- 
infrared SED is indeed that of an AO star where 
X^sSk^ttna^y^WHBitffl-ifeo^^ should find the 
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Fig. 2. — NGC 1514 as seen in various combinations of WISE and visible wavelength filters. The 
axes indicate arcminute offsets from the nominal central source position of a = 04 h 09 m 16!990, 
S = +30° 46' 33'.' 44 (J2000). a) WISE 3-color image, blue = W2 [4.6 pm], green = W3 [12 fim] and 
red = W4 d 

econ [22 /iinj, normalized to bring out the different infrared features, b) Visible 3-color 
image derived from POSS B, R, and I images, blue = B, green = B+R, red = R+I; rebinned 
to WISE's pixel scale and normalized so that the field stars are white, c) Visible wavelength 
combination of blue = POSS B, green = [O III], and red = Ha, normalized so that the nebulosity 
at (0.2, -0.6) is white, d) Combination of blue = POSS B, green = Ha, and red = W3, also 
normalized so that the nebulosity at (0.2, -0.6) is white. From c) and d), we conclude that the 
inner shell is remarkably "gray" across all wavelengths, in marked contrast to most other PNe. 
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Table 1: Flux of the Central Source of NGC 1514 

Observed Dereddened 
Band A Mag Extinc. Magn. Flux Magn. Flux 
Flux Corr. 





(jnn) 


(Jy) 


(mags) 


(mags) 


(Jy) 


(mags) 


(Jy) 


J 


1.235 


1594 


0.442 


8.190 


0.844 


7.748 


1.268 


H 


1.662 


1024 


0.275 


8.098 


0.590 


7.823 


0.761 


K 


2.159 


666.7 


0.133 


8.002 


0.420 


7.869 


0.475 


Wl 


3.353 


306.7 


0.058 


7.918 


0.209 


7.860 


0.220 


W2 


4.603 


170.7 


0.034 


7.907 


0.117 


7.873 


0.121 


W3 


11.561 


29.04 


0.048 


7.576 


0.027 


7.528 


0.028 


W4 


22.088 


8.239 


0.026 


> 7.5 


< 0.008 







same magnitude value at all wavelengths. Tak- 
ing the mean of the magnitudes from H through 
W2, the dispersion is reduced from 7.98 ± 0.08 
for the uncorrected case (where there is a ob- 
vious slope to the data) to 7.86 ± 0.02 after 
dereddening. The literature shows a fair spread 
in the reported B and V band values (though 
B — V is pretty consistently between 0.50 and 
0.55), but an average of a half dozen reported 
values shows that V = 9.53±0.09. Applying our 
value of Ay to this mean corrects the V magni- 
tude to 7.93, only 7% off the 7.86 found from the 
infrared magnitudes. In subsequent discussion, 
we assume that the best estimate of the intrinsic 
V magnitude of the CSPN is my = 7.86. 

3.2. Surface brightness 

Contour maps of NGC 1514 at the four WISE 
wavelengths are presented in Figure HI The neb- 
ular emission is clearly rising rapidly at longer 
wavelengths. The surface brightness of the neb- 
ula in the region of the bubbles roughly follows a 
power law F u oc i/ -2,0 for the four WISE bands. 
(This holds only for the WISE bands; the in- 
ner shell is brighter in the 2MASS bands than 
this power law would predict. Our assumption 
is that shorter wavelengths are dominated by 
line emission and perhaps scattered starlight.) 
It is also the case that the emission from the 
rings is rising even more rapidly (is redder) than 
that in the inner core, following a power law of 
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Fig. 3. — The SED of the central source. 
The plot shows a composite spectrum of an 
05V template (HD46223), an A0III template 
(HD123299), and a 9800 K blackbody, red- 
dened with Ay = 1.6 (dashed line). This is 
overlayed on the WISE and 2MASS data (open 



squares), anc 



a selection of UBV and near-IR 



photometry ( 


Arkhipova 1969. 


Kohoutek 


1967, 


Liller & Shao 




1968. iTaranova & Shenavrin 


200l 


Tvlenda et al.l 


1991). Except for the J 



flux, the SED of the CSPN is well represented 
by this composite spectrum and gives us 
confidence that the intrinsic visual magnitude 
of the central source is 7.86. 

Therefore, the nature of the nebular emission 
at mid-infrared wavelengths appears to be pri- 
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Fig. 4. — Surface brightness contour plots of 
NGC 1514 at the four WISE wavelengths. Con- 
tours for Wl and W2 range from to 0.8 (black) 
and 1.0 to 1.8 (white) MJy sr _1 in steps of 0.2; 
for W3: to 8 (black) and 10 to 18 (white) 
MJy sr _1 in steps of 2; and for W4: to 16 and 
20 to 32 MJy sr" 1 in steps of 4. 

marily thermal, even in regions where the line 
strength is relatively high. Though the emission 
from the brightest mid-infrared emission line, 
that of the [O IV] line at 25.89 /im, is reason- 
ably strong (see Section l3~4l) . it lies near the long 
wavelength cutoff of the W4 passband where the 
transmission is low, and thus cannot account for 
a substantial fraction of the observed flux. 

Figure [5] shows a plot of the total inte- 
grated flux of the nebula (all stars, including 
the CSPN, have been removed by profile fitting) 
from 1 to 100 /im, measured using a 5 arcminute 
diameter aperture for the 2MASS, WISE, and 
NVSS data, 3.5 arcminutes for the ISO/PHT 
data (the maximum extent available in the re- 
leased data), and the IRAS point source cata- 
log's values where we assume all of the flux will 
be captured in a single beam. These values are 
also listed in Table [2j 



At wavelengths shorter than 4.6 /im, the 
CSPN has a higher flux than the whole of the 
nebula. But at wavelengths longer than 4.6 /im, 
the nebula's flux density rises dramatically. A 
180 K blackbody curve has been drawn into the 
figure to guide the eye, but it is clear that mul- 
tiple temperature components are required to 
completely fit the data. A more detailed ther- 
mal model is not warranted given the few data 
points, but it is the case that the bulk of the 
emission comes from material that is quite cool, 
T < 200 K, and some perhaps as low as 30 K 
to account for the far- infrared fluxes. 
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Fig. 5. — A plot of the integrated nebula 
SED using the dereddened 2MASS, WISE, 
ISO/PHT, and IRAS fluxes. The central source 
flux (found by PSF fitting) is shown for compar- 
ison and is represented by the open circles, the 
integrated nebular flux (see text) is given by the 
filled circles, and the ISO and IRAS point source 
fluxes are indicated by the filled squares and tri- 
angles. An estimate of the contribution of the 
atomic lines was subtracted from the W3 and 
W4 fluxes, and this net thermal emission is in- 
dicated by the filled diamonds. A 9800 K black- 
body curve is overlayed on the CSPN fluxes 
while a 180 K blackbody is drawn through the 
WISE line-corrected W3 and W4 fluxes to sug- 
gest that the character of the mid-infrared neb- 
ular emission is thermal in nature. 
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Table 2: Infrared Integrated Nebula Fluxes 

Source/ Wavelength Observed Dereddened 
Band Flux Density Flux Density 





(pun) 


(Jy) 


(Jy) 


J 


1.24 


0.111 


0.166 


H 


1.66 


0.080 


0.103 


K 


2.16 


0.087 


0.098 


Wl 


3.35 


0.119 


0.126 


W2 


4.60 


0.148 


0.153 


W3 


11.56 


2.529 


2.644 


W4 


22.09 


8.264 


8.461 


ISO/PHT 


60 


6.205 




ISO/PHT 


90 


9.775 




IRAS 


60 


10.17 




IRAS 


100 


21.97 





3.3. Background 

There is a moderately bright, dusty back- 
ground presen t in th e region near NGC 1514; 
Schlegel et al. ( 19981 ) (as implemented by the 
NED Extinction Calculator) predict an Ay of 
2.252 magnitudes in this direction (/ = 165.53, 
b = -15.29). It appears that the PN is super- 
imposed on top of one of the brighter ridges of 
this background; we estimate a 12 //m bright- 
ness of 0.34 MJy sr _1 at the PN above the dark- 
est areas of the image. (The zero point is not 
yet calibrated, so we cannot know the absolute 
brightness.) A 3°x3° W3 field is shown in Fig- 
ure [6l The equivalent W4 field shows all the 
same structures visible at W3, while a compari- 
son field taken 45° farther off the galactic plane 
(I = 165.53, b = -60.29), where the predicted 
Ay is only 0.076, shows no such background 
structure. We are therefore confident that this 
background structure is real. 

The extended emission outside a 5 arcminute 
diameter circle around NGC 1514 is dominated 
by this striated background, and we see no ob- 
vious correlation between these structures and 
the nebula, either in position or symmetry. The 
PN is simply superimposed onto one of the 
denser portions of this region. This conclusion 
calls into question some of the results found 



by IWeinbergerl (jl999h and lArval et all lj20ld ) 
where they reported structures with degree an- 
gular scales that were linked with NGC 1514 
in the IRA S data at 12 fj,m. The emission 
Arval et al.l see in their Figure lc matches that 
of our Figure [6j but the much poorer spatial 
resolution of IRAS has smeared the background 
variations into the ring-like structure those au- 
thors report. Indeed, if we bin the WISE data 
down to two arcminute pixels, the appearance 
is similar enough to suppose that there is a 
link. However, the factor-of-60 higher resolution 
WISE data show that the background nebulos- 
ity is not correlated with the PN in any simple 
way. 

3.4. Spitzer spectra 

The Spitzer Space Telescope observed 
NGC 1514 only once as part of a wavelength cal- 
ibration program for the IRS instrument. Ac- 
cordingly, the exposures are very short, ranging 
from 15 seconds in the low-resolution data to 
60 seconds in the high resolution data, so that 
the overall signal-to-noise ratio is small. In ad- 
dition, the pointings between the four different 
spectrometer modules do not overlap, so that 
only on the central source do they all sample 
the same region. These data are publicly avail- 
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Fig. 6. — A 3x3 degree field in the vicinity of 
NGC 1514 shows the striated background visi- 
ble in the W3 band. The ridge upon which the 
PN is projected has a brightness of about 0.34 
MJy sr _1 above the darkest area of the region 
(which may not be zero flux). The PN itself is 
at the center. The white circles mask areas with 
known image defects: detector latents and opti- 
cal reflections, both caused by very bright stars. 
The equivalent W4 image, except for the factor- 
of-two lower spatial resolution, has all the same 
structures that are visible here. 

able and were accessed through the Spitzer Data 
Archive. 

Low resolution spectra were extracted for the 
central position by using the most central of 
the pointings while using the most off-center 
pointings to subtract the background spectrum. 
This spectrum is presented in Figure [71 Several 
facts are immediately apparent: 1) the short 
wavelength end of the IRS spectrum is domi- 
nated by emission from the 9800 K photosphere; 
2) the long wavelength end is dominated by 
100 K thermal emission from circumstellar ma- 
terial; 3) forbidden atomic lines from the nebula 
are quite strong even in just the line of sight 



to the CSPN; and 4) there is a complete lack 
of the PAH emission bands. Although PAH 
emission is usually dominant in nebulae, many 
PNe are known to have no PAH band emission 
and are instead dominated by H2 line emission; 
NGC 1514 appears to be a member of this class. 
Since the WISE W3 band encompasses all the 
PAH emission bands from 7.6 to 13.5 /xm, the 
Spitzer data suggests that all the structure seen 
in W3 arises from the same material that pro- 
duces the optical nebula, along with thermal 
emission from dust that is distributed through- 
out. 
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Fig. 7. — Spitzer low resolution spectrum at 
the position of the CSPN. Both stellar and neb- 
ular contributions are included in the flux. The 
dashed line represents a 9800 K blackbody pho- 
tosphere plus a 100 K dust continuum. The 
WISE relative spectral response curves (Wl is 
off scale to the left) are plotted below the spec- 
trum for reference; the scale is logarithmic from 
1% to 100%. Only the six labeled forbidden 
atomic lines are detected; the 10.51 /im [S IV] 
and 15.56 /im [Ne III] will contribute to the W3 
flux, but the 25.89 fim [O IV] line will contribute 
little to the W4 flux since the response is falling 
off rapidly. No significant PAH emission is seen. 



Only six emission lines are unambiguously 
detected at either spectral resolution: the 
25.89 /um [O IV] line is quite strong; the 
10.51 jum [S IV], 15.56 fan [Ne III], and 36.01 p 
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Fig. 8. — Spitzer high resolution spectra. Only 
these six lines are detected in the very short ex- 
posures. The spatial variations in intensity of 
five of the six lines match the intensity profile 
of the broadband W3 filter except at the center 
where the CSPN dominates the broadband pro- 
file (an arbitrarily-scaled version is plotted as 
the dotted curves) . Only in the [O IV] line does 
the distribution vary from the profile; it is con- 
centrated in the central 1-arcminute-diameter 
region. 

[Ne III] lines are fainter, but still well detected; 
and the 18.71 /im and 33.48 fim [S III] lines are 
very faint, but definitely present. The presence 
of these lines, in particular that of [O IV], in- 
dicates a high excitation nebula (central source 
of order 100,000 K), though not as high as 



some because of the lack of the [Ne V] line 
at 24.32 fj,m. 

Nine spatial positions along the PN's polar 
axis were sampled at high spectral resolution 
with the Short-High module, and seven posi- 
tions were obtained along the equatorial plane 
with the Long-High module (Figure [HJ). We note 
that, with the exception of the [O IV] line, all 
the lines vary in intensity as a function of posi- 
tion in the same manner as the broadband W3 
flux; the W3 intensity profile matches the en- 
velope of the line peaks surprisingly well. This 
implies that at least the sulphur (doubly and 
triply ionized) and neon are well mixed with the 
dust throughout the inner shell. The oxygen, 
on the other hand, (at least the triply ionized 
species) is concentrated in the inner 30" radius 
of the nebula, and is almost entirely absent out- 
side that. 

The outer two positions taken with Long- 
High module (75 arcsec from center) are the 
only two positions that would be dominated by 
flux from the rings. In both positions, the lines 
strengths (25.89, 33.48, and 36.01 /um) are much 
less than for the more central positions. Given 
that the nearest WISE band, W4, shows that 
the rings are brighter than the inner shell at 
long wavelengths, this suggests that lines are 
responsible for little of the ring emission; the 
rings must be almost entirely thermal emission. 

Because no pair of emission lines from the 
same species was observed at the same spatial 
position, we can only make crude estimates of 
line ratios. One of the better indicators of elec- 
tron density is the [S III] 18.7/33.5 /im line ra- 
tio. If we make the assumption that the un- 
derlying nebular brightness is roughly constant 
50" from the CSPN, then the ratio is a very 
low 0.33 (2.2 x 10~ 17 W m" 2 vs. 6.6 x 10~ 17 
W m~ 2 ). Even if our uniformity assumption is 
off by more than a factor of two either direc- 
tion, the inferred electron density is still well 
under 10 3 cm -3 . This agrees with the estimate 



by iKohoutekl (jl967l ) who derived a value of 290 



cm 3 for the inner shell and 140 cm 3 for the 
outer shell based on H/3 line strengths. 
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4. Discussion 

Given that the central source of NGC 1514 
is a binary system, the presence of these rings, 
while perhaps surprising in their form, is not en- 
tirely unexpected. Binarity and bipolar struc- 
tu res in PNe go h and-in-hand (e.g. the review 
bv lde Marcoll2009l . and references therein). The 
new infrared data, and particularly the discov- 
ery of the rings, allow us to refine our under- 
standing of NGC 1514. 

4.1. Distance 

The distance to NGC 1514 is quite uncer- 
tain. From the Hipparcos catalog of parallaxes, 
one can derive a distance of 185 ± 58 pc , while 
the st atistical methods summarized by IZhang 
(1995) can yield distances up to 1300 pc. Our 
WISE results do not shed a great deal of new 
light on the distance, but we can use our more 
confident estimate of my = 7.86 to help set 
plausible limits on the range of possibilities. 

In deriving the distance modulus to the neb- 
ula, the absolute magnitude of the CSPN has 
a larger uncertainty by far than the dereddened 
apparent visual magnitude. If the brighter com- 
ponent of the CSPN i s taken to be an AOIII gi- 
ant as determined bv lKohoutekl (119671) with an 
absol ute magnitude of ~ — 1 (jJaschek fe Gomez 
19981 ). then the inferred distance is ~ 600 pc. If, 
instead, the CSPN is a horizon tal branch A star 
as found by iGreensteinl (119721) with a n absolute 
magnitude ~ +1 (jde Boer et al.lll997l ). then the 
distance is ~ 240 pc. 

It seems likely, therefore, that distances 
> 600 pc s uch as those foun d by the H/3-to- 



radio ratio ( Cahn et al 



1992) or diameter-to- 



radio flux relationships (jDaubl Il982j ) are over- 
estimates. The Hipparcos distance (at +la) 
is broadly consistent with our lower luminosity 
estimate, thus we consider a distance of 200- 
300 pc quite plausible. 

4.2. Ring Properties 

For purposes of extracting geometric infor- 
mation, the rings can be modeled as two par- 



allel, unresolved rings (Figure [9j) with a depro- 
jected separation of 41" (0.05 pc, if we assume 
d = 250 pc), the southeast of which is 173" in 
diameter while the northwest is 177" in diam- 
eter (roughly 0.2 pc). The rings are tilted 59° 
from pole-on and rotated to a position angle of 
131°. A formal least squares fit using a downhill 
simplex method was attempted, but the irregu- 
larities in the ring brightness and the complex- 
ity of the interior shell morphology prevented a 
believable solution from being found. The val- 
ues quoted are estimated from "by eye" fits, but 
numerous trials showed noticeable errors in the 
diameter and separations at the ±1" level, in 
PA at ±0.5°, and in tilt at ±2°. Higher spa- 
tial resolution imaging than can be provided by 
WISE will be required to refine these estimates. 



W3 - smoothed Ha 
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Fig. 9. — The residuals of the ring model fit for 
the W3 image. A smoothed version of the Ha 
image has been subtracted to remove the inner 
shell from the W3 image (left panel) ; this image 
was then rotated clockwise 131° so the polar 
axis is vertical. The ring model (center) was 
then subtracted to reveal the residuals (right). 
While not a perfect fit due to the over-simplicity 
of the model, the general orientation and aspect 
ratios seem to be correct. 

By subtracting the ring models from the inte- 
grated flux measurements, we estimate the per- 
centage contribution of the rings to the total 
nebular flux. The rings contribute < 10% at Wl 
and W2, ~ 15% at W3, and > 30% at W4. The 
increasing fraction is consistent with the redder 
color of the rings noted in Section 13.21 and sug- 
gests a ring color temperature of ~ 160 K. 

This ring model fit can then be overlayed on 
all the images in order to see how the various 
features visible at differing wavelengths com- 
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pare (Figure [TO]) . Starting with the middle row, 
the fit falls well on the rings seen at all four 
WISE wavelengths (Wl is very similar to W2 
and has been omitted). At visible wavelengths, 
the ring fit falls entirely within the faint outer 
shell, but no evidence can be seen for the rings 
themselves at these wavelengths. In the far- 
infrared, the nebula is not well resolved, but it 
appears that it is elongated roughly along the 
equatorial plane of the model, suggesting that 
there is an extensive, dusty disk lying in the 
plane. This disk cannot be very dense given 
that the extinction to the central source is not 
large, and there are no color gradients across 
the nebula. 

At 1.4 GHz in the NVSS survey, the nebula 
is resolved into two spots correspondin g to the 
visible bubbles. Pazderska et al. ( 20091 ) report 
a total flux of 60 mJy at 30 GHz for NGC 1514. 
They note that this PN shows no sign of a 
high-frequency excess such as that attributed 



to em ission from spinning dust (jCasassus et al. 
20071 ). The radio emission from NGC 1514 is 



therefore purely free-free; this again suggests 
that the bubbles contain a significant quantity 
of ionized gas while the rings do not. 

4.3. A Physical Model, Refined 

Muthu fc Anandaraol fa00$ ) present a phys- 



ical model of NGC 1514 in their Section 4 that 
agrees very well with our observations. The 
overall ellipsoidal shape of the PN is still appar- 
ent in the infrared. The presence of the rings 
confirms the axisymmetric structure they as- 
sume. We refine the estimate of the position 
angle of the polar axis from 135° to 131°, and 
we establish a tilt angle of 31° out of the plane 
of the sky, agreeing with them that the south- 
east part of the nebula is the side that opens 
toward us. While they could not determine the 
relative angles between the bubble axis and the 
shell axis, we can confirm that they are coinci- 
dent (see, for example, Figure [TUj). This general 
configuration now seems quite secure. 

The formation of the rings themselves is still 
uncertain. At the simplest level, the rings are 



part of the dust that is associated with the 
photon-dominated region (PDR) that is, in 
turn, wrapped around the ionized zone. But, 
lacking specific kin ematic information on the 
rings (recall that iMuthu &: Anandaraol 12003 



concentrated on the inner shell, and where their 
sample points might have contained the rings, 
the rings are not detected) , we can only compare 
NGC 1514 to other objects that have similar 
morphologies. Two such objects are Hen 2-104 
and MyCn 18, although in both of these cases 
the structures are seen at visible wavelengths. 

The core of Hen 2-104 bears the most striking 
resemblance to NGC 1514, though Hen 2-104 it- 
self is a symbiotic star rather t han a PN and 
i t sho ws multiple sets of rings 



2001 



Corradi et al 



and updates bv lSantander-Garcia et al 



(2008), develop spatiokinematic models to 
explain the rings and accompanying struc- 
tures. They invok e the model proposed by 



Soker fc Rappaportl (|2000h that uses interact- 



ing winds from the central binary source to ex- 
plain the multiple rings and overall collimated 
structure. Such a model could partially explain 
NGC 1514 if we assume that the rings we ob- 
serve represent a particularly large mass loss 
event and are thus merely the brightest ones in 
the system, and that other rings and a more 
extensive structure may be detectable if greater 
sensitivity were available to us (and it were not 
washed out by the dusty background emission). 
However, the focusing mechanism of the winds 
cannot be so efficient as to prevent the forma- 
tion of the relatively uniform outer shell seen at 
visible wavelengths. 

The formation of MyCn 18 may follow similar 
lines. It has a very well defined hourglass shape 
with numerous pairs of parallel rings, though 
the overall structure is much more focused than 
either Hen 2-104 or NGC 1514. ISahai et al 



(1999 ) propose a mechanism related to that 
of ISoker k, Rappaport (|2000h . except that the 
rings are formed by a more collimated outflow 
jet interacting with a round circumstellar en- 
velope. This could again fit NGC 1514 if we 
assume that there was one event that was far 
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Fig. 10. — A comparison of a double ring model (drawn as the alternating black and white dotted 
curves) to the POSS B, O III, and Ha data (top row, left to right), the W2, W3, and W4 images 
(middle row), and the ISO 60 fj,m, ISO 90 fim, and 1.4 GHz NVSS maps (bottom row). It is clear 
that the rings are contained entirely within the outer shell. 
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stronger than any other, and any other events 
are simply too faint for WISE to detect. 

This still leaves the question of why the rings 
are visible only in the infrared. Our interpreta- 
tion is that, based in the spectral index of the 
ring emission discussed earlier, the rings are still 
present but are just too faint in comparison with 
the more uniformly distributed line emission in 
the outer shell. The rings are simply washed 
out at visible wavelengths. 

5. Conclusion 

The discovery of axisymmetric rings in 
NGC 1514 completes our perception of it as 
a planetary nebula formed from an aging bi- 
nary star system. No longer can NGC 1514 
be considered a simple planetary nebula "of a 
circular form" as Herschel perceived it. While 
the morphology is very complex, with numer- 
ous bubbles contained within the inner shell and 
rings contained within the outer shell, it joins 
the family of hourglass-shaped nebulae. 

It is likely that more such unexpected struc- 
tures will be found in PNe as the WISE survey 
data becomes widely available. As new hour- 
glass objects are found, the formation models 
can be put to more rigorous tests, and we will 
gain a deeper understanding of how these mys- 
terious and beautiful objects have come to exist. 
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